Weak Localization of Dirac Fermions in Graphene 
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In the presence of the charged impurities, we study the weak localization (WL) effect by evaluating 
the quantum interference correction (QIC) to the conductivity of Dirac fermions in graphene. With 
the inelastic scattering rate due to electron-electron interactions obtained from our previous work, 
we investigate the dependence of QIC on the carrier concentration, the temperature, the magnetic 
field and the size of the sample. It is found that WL is present in large size samples at finite carrier 
doping. Its strength becomes weakened/quenched when the sample size is less than a few microns 
at low temperatures as studied in the experiments. In the region close to zero doping, the system 
may become delocalized. The minimum conductivity at low temperature for experimental sample 
sizes is found to be close to the data. 
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Study of graphene is currently a focused area. The 
observations on the electronic transport properties show 
there is no weak localization (WL) of electrons in the ex- 
perimental samples P, [2| . This result seems to conflict 
with the general theory that the electrons should be lo- 
calized in two dimensional systems 0, ■ There have 
been theoretical works qualitatively discussing this issue 
primarily based on the impurity scatterings with zero- 
range potentials 0, 0, H, On the other hand, it has 
been shown that the charged impurities with screened 
Coulomb potentials [TlL IT^ are responsible for the 
observed carrier density dependence of the electric con- 
ductivity. Analysis of the WL effect in graphene using 
realistic model for impurity scatterings and its compari- 
son with experiments is still lacking. 

In this work, on the basis of the self-consistent Born 
approximation (SCBA) to the Dirac fermions under the 
charged impurity scatterings, we solve the integral ma- 
trix equations for the current-vertex corrections and the 
two-particle propagators. The WL is studied by evalu- 
ating the quantum interference correction (QIC) to the 
electric conductivity. With the inelastic scattering rate 
due to electron-electron interactions , the dependence 
of QIC on the carrier concentration, the temperature, the 
magnetic field and the size of the sample are investigated. 
The obtained results are compared with experiments. 

The Dirac fermionic nature of the electrons in graphene 
3.3) 3 ill have been observed by recent experiments 
ll.ll8|. Using the Pauh matrices ct's and r's to coordinate 
the electrons in the two sublattices of the honeycomb 
lattice and two valleys in the first Brillouin zone, respec- 
tively, and suppressing the spin indices for briefness, the 
Hamiltonian of the system is given by 

H = Y1 ^i^f" ■ ^^-^fc + E 4-,y^i1)i^k (1) 

k 



the valley indices, the momentum k is measured from 
the center of each valley, v {'^ 5.856 eVA) is the veloc- 
ity of electrons, V is the volume of system, and Vi{q) is 
the finite-range impurity potential. From the previous 
discussion for Vi{q) 12], we have 



VM = 



ni{-q)va{q)aa Jii{Q - q)viai 
ni{-Q ~ q)viai ni(--g)uo(g)cro 



(2) 



where ni{—q) is the Fourier component of the impurity 
density, vo{q) and vi are respectively the intravalley and 
intervalley impurity scattering potentials, and Q is a vec- 
tor from the center of valley 2 to that of the valley 1. 
Here, all the momenta are understood as vectors. The 
off-diagonal parts in Eq([2|) are slightly different from that 
in Refll2l where ao was used instead of txi because the ba- 



sis was vi = (cLi:cLi:cL2'4b2) with reflected y-axis 
in valley 2. But that basis is not convenient for dealing 
with the problem of WL. The potential vo{q) is given 
by the Thomas-Fermi (TF) type and vi ~ vo{Q) with 
Q — 47r/3a (a ~ 2.4 A as the lattice constant). It has 
been shown that the TF type potential gives rise to the 
same result for the electric conductivity of graphene as 
the one of the random-phase approximation [12] . The im- 
purity density chosen here is Ui — 1.15 x lO^^a^^, which 
fits the experimental data for the electric conductivity. 

Under the SCBA, the Green function G{k, uj) = [uj + 
fi ^ vk ■ — Y,{k,uj)]^^ and the self-energy Yi{k,uj) ~ 
TiQ{k, uj) + T.c{k, uj)k ■ (7Tz are determined self-consistently 
[3]. Here, ^ is the chemical potential, and k is the unit 
vector in k direction. The current vertex vTa{k,uji,uj2) 
(a — x,y) [Fig. 1(a)] can be expanded as 



Ta{k,UJl,UJ2) 



kq 



Y,yj{k,uji,uJ2)A'^{k) 

j=0 



(3) 



where = i4ai^4bi^4b2^4a2) is the electron opera- 
tor with a and b denoting the sublattice and 1 and 2 for 



where (A:) = r^aa, A"{k) = ctqit • k, (fc) = a ■ kaa, 
A"{k) = TzS-kuaS-k, andyj{k,uji,uj2) are determined by 
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FIG. 1: (Color online) (a) Current vertex, (b) Electric con- 
ductivity, (c) Quantum interference correction to the conduc- 
tivity, (d) Cooperon propagator. The solid line with arrow is 
the Green function. The dashed line is the impurity potential. 



four-coupled integral equations [I4I . The current-current 
correlation function [Fig. 1(b)] is obtained as 

P{uji,uj2) — ^yj{k,uji,uj2)Xj{k,uji,uj2) 
kj 

with Xj{k,oji,uj2) = TT[G{k,oji)A^{k)G{k,uj2)Af^{k)], 
for Lu's {oji and 1^2) = w ± iO = w^. 

To investigate the WL, we evaluate QIC [Fig. 1(c)] 
to the electric conductivity following the works of Refs. 
i and0. This correction is associated with two-particle 
propagator C^^^^^l^J^^^{k,k' ,q,u) (Cooperon). It obeys 
the Beth-Salpeter equation represented in Fig. 1(d). 
Here, the superscripts fs denote the valley indices, and 
the subscripts a's correspond to the sublattice indices. 
To solve this equation, we transform the Cooperons from 
the valley-sublattice space into the isospin-pscudospin 
space according to McCann et aZ.|9|, 



S1S2 » / , V^"si/'aia2 010203 



where Mj = EyE^AyA/, S's and A's are respectively the 
isospin and pseudospin operators defined as 



Ao = tqcto, 



Aa; — TxCTz, 



Sy — '^z^yi 

Ay — Ty(Jz , 



Az = r^o-Q. 



We will hereafter occasionally use the indices 0,1,2,3 for 
0,x,y,z, respectively. The equation for C^^g^^ reads, 

(fc, k',q) - (k, fcl , q) [Wi^ s' ( I fci - k' \)Sw 



+ Ci[Aki,k',q)] 
where 11^^^ is the element of matrix ft' defined as 
Il\k,ki,q) = W'{\k - ki\)h{ki,q) 



(4) 



(5) 



with W'^^,i\k - fell) = n,[v^{\k - fell) + vfiSio - 
Siz)i—^)^]Sss' , and the element of h is given by hsgi{ki, q) 
= Tr[G(-fc+,w+)EsG(-^r,w-)Sj,]/4 with kf = ki ± 
q/2. W^{\k — fc'l) is the isospin-pseudospin space rep- 
resentation of the interaction n^^od^ ~ ^'l)'^aia3'^a2'Q4 

ing the conjugate valley of j) shown as the dashed line 
with a cross in Fig. 1(d). The second term of this inter- 
action stems from the fact that when a particle is scat- 
tered from valley 1 to valley 2 another particle should be 
scattered inversely so that the total momentum of the 
Cooperon is unchanged. From Eq. ([3]), it is seen that 
the pseudospin of the Cooperon is unchanged during the 
impurity scatterings. We hereafter denote C"^, as C[g, . 

To get the solution to Eq. (jl]) , we firstly need to solve 
the eigenvalue problem of II^.^^ (fc, ki,q): 

^ ^ n[,^{k,h,q)^sAkl) - X\q)^s{k). (6) 



V ^ 

fcl, Si 



For I — and <f = 0, a solution is A"(0) — 1, and 

**(fc) = [Ao(fc,w),-Ac(fc,w)cos(/), -Ac(fc,w)sin(/i,0] 

where Ao(fc, lo) — ImEo(fc, ), Ac(fc, lo) = ImEc(/i;, i^^), 
and (/> is the angle of k. The four components of ^'(fc) 
correspond to s = 0, y, z respectively. The solution of 
A°(0) = 1 is the most important one which gives rise to 
the diverging contribution to the Cooperon. ^(fc) is a 
result of the SCBA for finite-range impurity scatterings. 
For the zero-range potential, only the first component of 
^(fc) survives and is a constant. One then needs to solve 
a scalar equation instead of the matrix equation. 

For g ^ 0, we consider only the case of small q since 
where QIC is significant. At small g, by expanding 
n°(fc, to second order in q and using perturbation 
treatment, we obtain 



\\q) 



1 



^ *t(fe)ftO(^^ g)«'(A;') = 1 - doQ^ 



kk' 



to the first order in the perturbation, where (^'I^E') = 
X^^, ^''''(fc)\l/(fc)/F and do is a positive constant. To the 
0th order, the eigenfunction is unchanged. 

Note that the difference between ft' for I ^ and ft" 
comes from the intervalley scattering term in . By 
regarding to this difference as a perturbation, we obtain 



\\q)^\\Q)-diq^ 



for 







to the first order, and the same eigenfunction as \E'(fc) to 
the 0th order in the perturbation. Here A' (0) is less than 
1, and approaches to 1 as the carrier density becomes 
extremely low. 

We now describe the solution C\^, to Eq. (g]). Sup- 
pose all the eigenfunctions of Eq. ([6]) have been ob- 
tained as ^'Ji(A:, q) (being a column vector in the isospin 
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FIG. 2: (Color online) Quantum interference correction Sg as 
function of electron doping concentration S. The inset shows 
the calculated (green circles) and the experimental (hollow 
circles) results for the electric conductivity. 



FIG. 3: (Color online) Quantum interference correction Sg as 
function of temperature T at 5 = 1.0 x 10"'', and 1.0 x 10~^ 
for sizes L — 1 fim and L > 10 /im. 



space) with eigenvalue X[^{q), n — 1, 2,- • • . Then, as 
a matrix in the isospin space, C' can be expanded as 

C'(fc,fc',g) = E„cUg)*Ufc,'7)*5!(fc',<7l. FromEq. Q, 
it can be shown that c5j(g) oc [1 — X^iq)] ^- There- 
fore, the predominant contribution to comes from the 
state with the lowest |1 — A'|. We will therefore take 
into account only the state of the lowest |1 — A'| for 
each I. The state for each / obtained above is just the 
one. We hereafter denote cjj(g) of the lowest state as 
c\q) EE q/[1 - \\q)] with Q = (^'|n'W^'|*)/((*|*))2. 

With the Cooperon C', the QIC to the current-current 
correlation function P(ti)^,u;+) is given by [Fig. 1(c)] 



5P{u;-,uj+) ^ ^Y,Tv[Z\k,q,Lo)C\~k,k,q)] (7) 

kql 



where Z {k,q,uj) is a matrix with elements Z^^,'s^ 

Zl, ^ Tr[V'{k+ ,iu+ ,uj')MlV{-k- ,uj- ,Lu+)Mlr], 

and V{k,LL)i,u!2) = G(k,uji)T rc{k,uji,uj2)G{k,uj2)- Here 
V*{k, uji, UJ2) is the transpose of V{k, lui, ^2)- By noticing 

Ml = M°Ai, AiM°A^ = for l = x,y,z 

and the operator A; commutes with G and F^,, we have 
Z^^i = —Z'^s' ^'^^ I — ^iVi ^- This result means that the 
QIC by the pseudospin singlet {I — 0) is negative, while it 
is positive by the pseudospin triplets {I ~ x,y,z). Since 
the g-integral in Eq. ([7]) comes predominantly from the 
small q region, the g-dependence of can thereby be 
neglected. Carrying out the g-integral, we get 

)-E"'^ln ^_,.(0) + ,,,2 . (8) 

'f''{k)Z°{k,0,uj)^{-k), (9) 



SttV ^ 

k 



where uq — —1, ni^x.y,z ~ 1, qo and qi are the lower and 
upper cutoffs of the (f-integral. 

The lower cutoff qo is given by qo — max(L,^^, L^^) 
where Lin is the length the electrons diffuse within an in- 
elastic collision time Tin and L the length scale of the sys- 
tem. We have recently studied the interacting electrons 
in graphene using renormalized-ring-diagram approxima- 
tion [13|. At very low doping {S < 1.0 x 10~^, the doped 
electrons per site), from the result for the self-energy, 
due to the inter-electronic Coulomb interaction is esti- 
mated as Tin ~ 0A62v/aT^, where T is the temperature. 
Lin is then given by Lin = ip^T^inl'^^^'^ 4] where the 
elastic collision time r is determined by the non QIC- 
corrected conductivity goj 7" = TT-T^go/vkpe^ (with kp as 
the Fermi wavenumber). For low carrier density, we find 
that Lin is about a few microns for 4 K < T < 20 K. On 
the other hand, the upper limit is qi — Lq with Lo = vt 
as the length of mean free path. 

In Fig. 2, QIC to the conductivity Sg « SP{Q- ,0+)/2^T 
a,t T = 10 K is shown as function of the electron doping 
concentration S, and compared with the corrected con- 
ductivity g obtained from Fig. 1(b) -I- Fig. 1(c) (for 
L = 1/im) and the experimental result [l|. For smaller 
size system L < Lin, Sg is larger with S fixed. This is 
because the lower cutoff go moves up for small L so that 
the negative contribution from the I = channel is weak- 
ened. Moreover, 6g increases with decreasing doping. 6g 
is positive at low doping. This stems from the fact that 
at small S the screening is weak, the Fermi circle and the 
typical momentum transfer q (~ 2kF) are small, leading 
to stronger vo{q) than vi. For weak vi/vo{q), all A'(0)'s 
(Z 7^ 0) close to 1, the contribution to Sg from each pseu- 
dospin channel has almost the same magnitude. After 
one of Z ^ is canceled by the I = channel, the net 
contribution to Sg is positive. Physically, Dirac fermions 
cannot be scattered to exactly the backwards direction 
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FIG. 4: (Color online) Magneto-resistance 5R{B)/R{0) as 
function of i? in unit of Tesla at T = 4 K for various sample 
sizes and electron doping concentrations. 



in case oi Vi = and the WL is absent. Therefore, the 
system may become delocalized at very low carrier con- 
centration. On the other hand, vifith increasing 6, the 
strength of the intervalley scatterings becomes stronger, 
leading to the appearance of WL in large size samples. 

In Fig. 3, we exhibit the temperature and size de- 
pendences of Sg. For system of large scale L > Lin, 5g 
shows the logarithmic behavior at low temperature. Such 
a temperature dependence is absent for small size system. 
With L < 1 /im scale of the samples, the WL effect could 
hardly be observed in experiment at low temperature [l[ . 

For the minimum electric conductivity the non 
corrected value is about 1.7 (with unit e^//i) [T^- At 
T = 10 K, the contribution of QIC is Sg^ = 3.3 for L = 
1 /.tm, and its magnitude decreases and saturates to 2 for 
L > 10/im as shown in Fig. 2. Therefore, gm is around 
4. On the other hand, the QIC varies with T as shown 
in Fig. 3. For L = 5/im, we find Sgm varies in the range 
1 < 6grn < 2.8 for 4 K < T < 20 K. As a result, for 1 
/im < L < 5/im (the experimental sample sizes) and 4 K 
< T < 20 K, is close to the experimental data [ij. 

In the presence of an external magnetic field B, the 
Cooperons are in the quantized Landau states [l3| ■ The 
(f-summation in Eq. ([7]) is then replaced with the sum- 
mation over those Landau levels (each of them with a 
degeneracy BV/tt). By using the same g-cutoffs as in 
obtaining Eq. dH), only those states of energy levels in 



the range {q^/'i., (7i/2) (in units of a 



1) need to be 



summed up. In Fig. 4, we show the relative magneto- 
resistance 6R{B)/R{0) as function of i? at T = 4 K and 
6 = 2.5 X 10^*^ for various sample sizes. The result for 
6 = 1 X lO^"' and L = 1 /im is also depicted for compar- 
ison. As seen from Fig. 4, the magneto-resistance very 
sensitively depends on the sample size. When L < Lin, 
the number of the Landau levels is determined by L. 
Since smaller sample confines less Landau states, the 
magnetic field effect is largely reduced in small samples. 



At very low 5, though the QIC to the electric conductiv- 
ity for B = appears like a delocalization effect as given 
in Fig. 2, the effect of magnetic field presents a suppres- 
sion of WL giving rise to reduced resistance as shown in 
Fig. 4. Such a magnetic effect is qualitatively consistent 
with the experimental observation 0]. 

In summary, using SCBA to the Dirac fermions under 
the charged impurity model, we have studied the WL 
effect in graphene by evaluating QIC to the electric con- 
ductivity. Since the model justified by the existing the- 
ories [lol 



111 . |12| is more realistic than the zero-range 



potential model, the present calculation should give rise 
to the reasonable predictions for WL in graphene. With 
the inelastic scattering rate [11], we have investigated 
the dependence of QIC on the carrier concentration, the 
temperature, the magnetic field and the size of the sam- 
ple. It is found that WL is present in large size samples 
at finite carrier doping. The strength of WL becomes 
weakened/quenched when the sample size < Lin (about 
a few microns at low temperatures) as studied in the 
experiment [l|, Close to region of zero doping, the 
system may be delocalized. In addition, the minimum 
conductivity obtained for experimental sample sizes at 
low temperature is close to the data. 
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